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Abstract 

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a 
validated, three-dimensional, large-eddy simulation model. Numerical simulations are performed in order to 
isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. 
Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an 
approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases 
clearly with increasing ambient turbulence level, which is consistent with field observations. 

Based on the results from the numerical simulations, simple decay models are proposed as functions of 
dimensionless ambient turbulence intensity ( 77 ) and dimensionless time (T) for the circulation averaged over 
a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex 
decay model is proposed for weak turbulence, while an exponential type of vortex decay model can be 
applied for strong turbulence, A relationship for the vortex descent based on above vortex decay model is 
also proposed. Although the proposed models are based on simulations assuming neutral stratification, the 
model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable 
atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in 
reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely 
underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. 
The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an 
important influence on vortex decay when ambient levels of turbulence are weak. 


1. Introduction 

For the purpose of increasing airport capacity, a 
system is being developed under the National Aero- 
nautic and Space Administration's (NASA) Termi- 
nal Area Productivity (TAP) program that will con- 
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trol aircraft spacing within the narrow approach cor- 
ridors of airports. The system, called the Aircraft 
Vortex Spacing System 1 ’ 2 (AVOSS), will determine 
safe operating spacing between arriving and depart- 
ing aircraft as based on the observed or predicted 
weather conditions. This system will provide a safe 
reduction in separation of aircraft compared to the 
now-existing flight rules, which are based on aircraft 
weight categories. In order to develop this system, 
research is being focused on understanding how air- 
craft wake vortices interact with the atmosphere. 
Previous studies indicate that transport and decay 
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of wake vortices are strongly affected by ambient at- 
mospheric parameters such as wind shear, stratifi- 
cation, and turbulence, as well as by proximity of 
aircraft to the ground. 3 In this study we focus on 
the effects of three-dimensional ambient atmospheric 
turbulence on vortex decay and descent using a val- 
idated large eddy simulation (LES) model. 

The mechanisms and causes for the decay of the 
wake vortices are somewhat controversial. 4 Some re- 
searchers believe that vortices generally do not decay 
until three-dimensional instabilities such as Crow 
instability 5 lead to their sudden destruction. 6 How- 
ever, this view is held in spite of the overwhelming 
observational evidence that vortices decay at differ- 
ent rates depending upon the ambient meteorology, 3 
and atmospheric turbulence plays a key role. 7,8 Us- 
ing a second-order turbulence closure model, Bilanin 
et al. 9 shows that the rate of decay of a vortex 
pair increases with increasing background dissipa- 
tion rate. Recent analyses of field observation data 
by Sarpkaya 10 clearly show that vortex circulation 
decays much faster in the presence of strong turbu- 
lence than in weak turbulence. The analyses of LES 
data 3,11,12 also reveal that three-dimensional atmo- 
spheric turbulence enhances circulation decay. On 
the other hand, the descent speed of a vortex pair 
due to mutual induction decreases with increasing 
circulation decay as long as the separation distance 
of a vortex pair remains constant. From data ob- 
tained in the water tank experiments, Sarpkaya and 
Daly 13 find that the descent speed of vortices de- 
creases with increasing dissipation rate of the ambi- 
ent turbulence (i.e., increasing ambient turbulence 
level), which in turn implies that stronger turbu- 
lence yields larger reduction in circulation. Stable 
stratification can also hasten vortex decay through 
baroclinic generation of counter-sign vorticity along 
the vortex oval’s periphery. 3,14,15 Since most of field 
observation data are subjected to the influence of 
both turbulence and stratification, it is difficult to 
distinguish their effects on circulation decay and vor- 
tex descent. In the present study, we consider only 
a neutrally stratified atmosphere in order to isolate 
atmospheric turbulence effects on vortex decay. In- 
vestigations of wake vortex decay in a non-neutral 
environment are ongoing and will be reported in a 
future paper. 

The enhancement of vortex decay due to am- 
bient turbulence has been an important factor in 
wake-prediction models. 7,8 Donaldson and Bilanin 7 
proposed a formula describing the relation between 


atmospheric turbulence and vortex circulation: 


dT 

dt 



(i) 


where T is circulation, q = ( u 2 + v 2 + w 2 ) 1 ^ is the 
turbulent velocity scale, and bo initial separation dis- 
tance of a vortex pair. Greene 8 incorporated Eq.(l) 
into his analytic model for the contribution of at- 
mospheric turbulence effects on wake vortex motion 
and decay. Equation (1) implies that vortex circula- 
tion decays exponentially and depends only on ambi- 
ent turbulence, represented by q which is related to 
the turbulence kinetic energy ( TKE = q 2 / 2). The 
TKE or q is primarily determined by large energy- 
containing eddies. According to Eq.( 1), then, large 
scale eddies are related to the decay of vortex circu- 
lation. Although, large energy-containing eddies can 
lead to a large deformation and transport of vortices 
as well as to the development of the large scale Crow 
instability, 3,16,17,18 it is the small scale (order of 6 0 ) 
turbulent eddies that are directly associated with the 
decay of vortex circulation. These smaller scales of 
turbulent eddies are likely to be more homogeneous 
and isotropic and mostly lie within the Kolmogorov 
inertial subrange which is characterized only by tur- 
bulent kinetic energy dissipation rate (e). 19 There- 
fore, one can use e rather than q as the turbulence 
parameter for controlling the vortex decay. The two 
are related as e = q 3 /l t , where l t is a characteristic 
large-eddy length scale, called the dissipation length. 

In fact, e has been commonly used in the 
literature 13,17,18,20,21 as an appropriate parameter to 
investigate the effects of ambient turbulence on the 
lifespan of wake vortices. The strength of turbulence 
is often represented by the dimensionless turbulence 
intensity q, defined as 


q = 


(<M 1/3 


( 2 ) 


where 


V 0 = 


( 3 ) 


27T&0 

Here, represents circulation at very large radial 
distance. Note that q represents the ratio of the 
characteristic turbulent velocity scale at the scale of 


the vortex separation distance ((e&o) 1 ^ 3 ) to the speed 
of descent of the vortex pair by mutual induction 

m 

Our studies are focused on the far-field, post roll- 
up stage of the vortex pair, in which each vortex pos- 
sesses a well-developed structure. For later discus- 
sion, we define nondimensional parameters for time 
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( t ), descent distance of the vortices (h), and radial 
distance (r) as 


T 


V 0 t 
bo ’ 


h r 

= and R= 
bo bo 


(4) 


In Section 2, we describe the LES model and the 
modifications required for initial conditions. In Sec- 
tion 3, we present the results from systematic numer- 
ical experiments in terms of nondimensional turbu- 
lence intensity and also propose new vortex models 
for vortex decay and descent, based on the LES re- 
sults. In Section 4, we apply the proposed models to 
available observation data. Finally, in Section 5, we 
summarize our results and draw some conclusions. 


2. The model and initial conditions 

2.1 The model 

The numerical model used in the present study 
is a three-dimensional, nonlinear, compressible, non- 
hydrostatic LES model, called the Terminal Area 
Simulation System 22 (TASS), which has been 
adopted for simulation of interaction of wake vor- 
tices with the atmosphere. 23,24,25 Recently, this code 
has been used to study the development of the Crow 
instability in an artificially generated homogeneous 
atmospheric turbulence. 18 The agreement between 
the calculated and measured maximum instability 
wavelengths and vortex lifespans provided a good 
validation of the present TASS LES code. The model 
equations are basically identical to those used in 
Han et al. 18 except for the modified subgrid model 
which accounts for the flow rotation effects. 

Turbulence is strongly affected by the rotation of 
the swirling flow of vortex. According to Rayleigh’s 
well-known stability criterion, 26 perturbations are 
suppressed in an axisymmetric vortex if the circu- 
lation is increasing with radial distance, but per- 
turbations would become unstable if the circulation 
is decreasing with radial distance. Since the circu- 
lation increases rapidly with radial distance in the 
vortex core region, any pre-existing turbulence is 
suppressed there and the vortex core^ expands very 
little with time. 15 In the present numerical model, 
a scheme for the subgrid eddy viscosity which takes 
account of the flow rotation effects is used to avoid 
unrealistic core growth. 


*The vortex core is defined as the radius of maximum tan- 
gential velocity. 


For subgrid turbulence, TASS currently uses a 
conventional Smagorinsky model 27 with modifica- 
tions for stratification effects as: 


K = (c s A) 2 \D\ (1 - aRis) 0 ' 5 , 

Ri s = - and TV 2 = | (5) 


where K is the subgrid eddy viscosity, Ris is the 
Richardson number due to stratification, A is the 
filter width, D is the rate of deformation, N is the 
Brunt- Vaisala frequency, g is the gravity, 0 is the 
potential temperature, and c s = 0.16 and a — 3 are 
constants. After initial sensitivity tests, we have re- 
alized that the Smagorinsky subgrid closure model 
overpredicts the decay rate in the vortex core re- 
gion and causes unrealistic growth of the core size. 
Since vortices have been observed to maintain their 
core strength and size for a much longer time, 10 it is 
necessary to incorporate the rotation effects into the 
Smagorinsky model. A simple modification of the 
eddy viscosity to include rotation effect analogous 
to the stratification effect has been developed. 28 

The modified subgrid eddy viscosity is given as 

K = (c a Af\D\{l-aRi s -(3Ri R f 5 . (6) 


Here Rift is the rotational Richardson number, de- 
fined as 28 

n 2 \n\ 

Rir ~ d 2 + \d\ ’ ^ 


where Q is the magnitude of three-dimensional vor- 
ticity. Equations (6) and (7) indicate that the eddy 
viscosity, AT, can be effectively reduced within the 
core of vortices due to large |0|/|D|. Since the above 
formula cannot discriminate between shear flow and 
flow with coherent rotation, a discriminator function 
is applied that turns the formula off (defaulting to 
the original Smagorinsky model) in the absence of 
coherent rotation. Details of this formula can be 
found in Ref. 28. 

Periodic boundary conditions are imposed at all 
domain boundaries. Domain size used in our simu- 
lations is (L x , L y , L z ) = (2.5&o, 5&o> 5&o)=(80Ax, 
160 Ay, 160Az) where b 0 = 32 m, grid size Ax = Ay 
= A z = lm, and x, y and z correspond to the axial, 
lateral and vertical directions of the vortex system 
and corresponding velocity components are u, v and 
w , respectively. The smaller domain size in the ax- 
ial direction (2.56o), which can save much computing 
time, suppresses the development of Crow instability 
of which theoretical maximum wavelength is about 
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Figure 1 : (a) Vertical velocity field and (b) its one-dimensional energy spectrum in a steady state turbulence 
before vortex injection. 


8.6bo, 5 and thus statistically homogeneous decay be- 
havior is anticipated along the axial direction. The 
domain size of 5&o in lateral and vertical directions 
is sufficiently large to minimize boundary influences. 
Due to limited computing resources and in order to 
allow the core to be resolved, the initial core size 
(r c = 4 m) is somewhat larger than the typical value 
observed behind aircraft (about 0.05 B 23 ). Ground 
effects axe not included in the simulations, thus re- 
sults are applicable only to out-of-ground-effect wake 
vortex behavior. 

2.2 Initial conditions 

Since we want to study the effect of ambient tur- 
bulence on the vortex decay and descent, it is of cru- 
cial importance to obtain an initially homogeneous 
and isotropic turbulence field. Toward this purpose, 
the initial turbulence field is allowed to develop un- 
der an artificial external forcing at low wavenumbers 
using a three-dimensional Fast Fourier Transform, as 
described in Han et al . 18 

Figure 1 shows the vertical velocity field and its 
one-dimensional longitudinal and transverse spec- 
trum with —5/3 slope of Kolmogorov’s spectrum 
when the turbulent flow field reaches a statistically 
steady state in which the mean turbulence kinetic 
energy fluctuates in time around a constant value 
due to viscous dissipation against external forcing. 
The integral length scale is calculated as L 33 « 11.4m 
and the large-eddy turnover time, defined as t e = 
^ 33 /('^ 2 ) 1 ^ 2 , is estimated to be about 51.6s. Here 


subscript 3 denotes vertical direction and ( ) repre- 
sents the domain average value. The isotropy pa- 
rameter /, defined as I\ = [(u 2 )/(w 2 )] 1 ^ 2 or I 2 = 

[{v 2 )/{w 2 )] , fluctuates only a few percent around 

its expected value of one for isotropic turbulence. 
Therefore, our simulated turbulence is close to sta- 
tistical isotropy. The TKE dissipation rate e is es- 
timated from the well known form of the spectrum 
in the inertial subrange and fitting the same to the 
simulated spectra in Fig. 1 . 

The initial vortex system is representative of the 
post roll-up, wake- vortex velocity field. A vortex 
model recently developed by Proctor 3 is adopted in 
our experiments. This model is empirical as it is 
based on field observations of several wake vortices 
measured early in their evolution. Its tangential ve- 
locity, V, is represented as 


V(r) = ^ (i-e- 10 ( r / B >°- 75 ), (8) 

where Too is the circulation at r >> r c ( r c is the 
initial vortex core radius defined as radial distance 
of peak tangential velocity) and B is the span of the 
generating aircraft. The values assumed for initial 
vortex separation and circulation are derived from 
aircraft weight (W) } wingspan (J3), air density (p), 
and airspeed ( V a ) according to the conventional as- 
sumption of elliptically-loaded wing, i.e., 



and 


Too = 


4 W 
irBpVa 


( 9 ) 
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Note that the velocity field in Eq.( 8 ) depends on the 
wingspan B instead of the core radius r c which may 
not be easy to measure with a good accuracy; Eq.(8 ) 
is applied only at r > r c . For r < r c , the model is 
matched with the Lamb model , 29 i.e., 

v<r) = 

^_ e -1.2527(r/r c )^ {10 ) 

Normalized circulation, tangential velocity, and vor- 
ticity from the above model are shown in Fig. 2 as 
function of the normalized radial distance R = r/ 6 0 , 
assuming r c /bo = 0.125. 



Figure 2: Initial circulation, tangential velocity, and 
vorticity for each vortex normalized by their maxi- 
mum values. 

A two-dimensional counter-rotating vortex pair 
prescribed by Eq.(8) and Eq. (10) is initialized uni- 
formly along the axial direction at the time when 
the ambient turbulence reaches a steady state after 
integration. In order to maintain the vortices in the 
middle of the domain, the grid is allowed to translate 
downward at the speed of 90% of Vo- Nondimen- 
sional ambient turbulence strength, 77 , is obtained 
by varying the circulation rather than by varying e 
to save computing time and thus all the following 
results are presented as nondimensional forms. As- 
suming that Too = 400m 2 s _1 and bo = 40 m, which 
are typical for DC — 10 aircraft, the range of typ- 
ical values of 77 in the lower atmosphere is about 
0 . 01 - 0 . 5 , with the maximum in the early afternoon 
and the minimum after midnight 3 . Simulations are 


conducted for six values of 77 , which are within this 
typical range. These cases are divided into three 
turbulence strength groups: 77 = 0.0302, 0.0702 for 
weak turbulence, 77 = 0.1006, 0.1509 for moder- 
ate turbulence, and 77 = 0.3018, 0.5031 for strong 
turbulence, where “weak” represents the strength 
of typical boundary layer turbulence during stable 
night-time conditions and “strong” represents typi- 
cal strength during sunny daytime conditions. In or- 
der to isolate atmospheric turbulence effect on vor- 
tex decay, the numerical simulations presented in 
this paper assume neutral stratification. 

3. Numerical results 

External forcing for maintaining ambient turbu- 
lence strength is deactivated to save computing time 
when the initial vortices are injected, but the differ- 
ence in the result for activating or deactivating the 
forcing has been found to be negligible (not shown). 
This confirms the arguments by Corjon et al . 30 that 
the time scale of the ambient turbulence compared 
to that of the vortex is sufficiently large to obtain the 
main characteristics of the effects of ambient turbu- 
lence on the wake vortices. 

During simulations, the vortices descend smoothly 
without any development of large scale instabilities 
such as Crow instability, which is not allowed due to 
the small domain size used in the axial direction. 

The circulation at any radial distance r can be 
easily determined from the area integral of the ax- 
ial vorticity (C) over a region defined by the radial 
distance, i.e., 

r r = JJ<dA. (11) 

Obtained at each radial distance are a total 160 cir- 
culation values for a pair of vortices in 80 y — z planes 
from which we have calculated averages and stan- 
dard deviations. As shown in Table 1 , the circula- 
tion fluctuations due to ambient turbulence increase 
with increasing turbulence level. The standard de- 
viations appear to be small for most of turbulence 
levels (less than 10 %) at least at the end of the 
simulation time, except for the strongest turbulence 
strength of 77 = 0.5031 for which the values are more 
than 20 %. In the following, we consider only the av- 
erage of circulations to deduce an appropriate model 
for vortex decay and descent. 

3.1 Circulation decay 
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1 

T 

3r c 

4r c 

5 r c 

0.0302 

7 

2.8 

2.0 

2.1 

0.0702 

6 

3.8 

4.1 

4.3 

0.1006 

6 

3.5 

4.5 

4.6 

0.1509 

6 

6.6 

9.0 

9.9 

0.3018 

3 

10.1 

9.6 

8.5 

0.5031 

3 

31.4 

28.2 

23.1 


Table 1: Standard deviation of the circulation rela- 
tive to its mean value (%) at varying radial distances 
and dimensionless turbulence strength 77 at the final 
dimensionless simulation time (T). 


Figure 3 shows the decay of circulation with time 
for varying radial distances as well as varying ambi- 
ent turbulence levels. As shown in Fig. 3, the decay 
rate of the circulation increases clearly with increas- 
ing ambient turbulence level and appear to decrease 
with increasing radial distance, which is consistent 
with field observations. 10 

In order to investigate the vortex decay behav- 
ior in more detail for varying radial distances, the 
circulation evolution with a best-fitted function for 
each radial distance is plotted in Fig. 4 for relatively 
weak turbulence (77 = 0.0702), in Fig. 5 for relatively 
moderate turbulence (77 = 0.1509), and in Fig. 6 
for strong turbulence (77 = 0.5031). As evidenced 
in Figs. 4-6, the decay of the circulation appears 
to follow a Gaussian function, e“ aT , for weak and 
moderate turbulence (hereafter model G), while for 
strong turbulence, it appears to follow an exponen- 
tial function, e -6T , at smaller radial distances (here- 
after model E) but a Gaussian function at larger 
radial distances, where a and b would be functions 
of 77 and radial distance. Although partly shown in 
Figs. 4-6, the circulation at radial distances larger 
than 0.6 b 0 appears to follow the Gaussian decay for 
all levels of turbulence. An exponential decay for- 
mula, Eq.( 1), was first proposed by Donaldson and 
Bilanin, 7 with q rather than e as the turbulence pa- 
rameter controlling the vortex decay. The Gaussian 
type of the decay (model G) replicates the much 
slower decay during early evolution of the vortices. 
This slow decay at earlier time has also been re- 
ported in Corjon and Darracq’s 31 LES results, in 
which the circulation begins to decay after a period 
of about 10 s. Based on our LES results described 
above, two types of models (i.e., model E and model 
G) are proposed for vortex decay and their formula- 
tions are given in Appendix A. 

Previous investigators have quantified the 


strength of a vortex in terms of a parameter called 
average circulation, which is defined as 23,32 

nr dr 

r n ,r 2 = ( 12 ) 

Jri 

where 7*1 and r<x are the radial distances of the aver- 
aging interval. This parameter is desirable, since 
it relates to the rolling moment of an encounter- 
ing aircraft and provides a more stable measurement 
than the local circulation . 23,32 In particular, the av- 
erage circulation within 15m from the vortex cen- 
ter has been considered to be closely related to the 
hazard . 33,34 In the present study, we compute the 
average circulation from r = 0.4 60 to r = 0.6 6 o- 
The lower limit is large enough in order to avoid 
any artificial contamination near the core. The up- 
per limit is small enough to minimize the effect of 
the accompanying vortex. 

Since T is a function of R , it is difficult to obtain 
decay formula for the average circulation directly 
from Eq.(AlA) and Eq.{A19) based on the definition 
of the average circulation, i.e., Eq. (12). Instead, for 
the decay of the average circulation, Eq.(AlA) and 
Eq.(A19) are approximated as 


dT 

dT 


turbulence 



for model E (13) 


and 


df 

dT 


turbulence 



for model G, (14) 


where c\ and C 2 are empirical constant and R is a 
value at the middle point of an average range of ra- 
dial distances (e.g., R becomes 0.5 for the average 
range of 0.46 o -0.66 o used in the present study). Fur- 
thermore, integrations of £g.(13) and Eq.(14) with 
respect to dimensionless timell (assuming ambient 
turbulence is only the process affecting vortex de- 
cay) yields 


f = f 0 e-( Cl " / " 2)T for model E (15) 

and 

f = f 0 e ~( C2rj2 / R2 } T 2 for model G. (16) 

As can be seen from Fig. 7, the model predictions 
of average circulation agree much better with the 
LES data than those of local circulations, showing 

^It is assumed that 77 does not change with time. 
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Figure 3: Decay of circulation with time for varying radial distances and ambient turbulence strengths. 
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Figure 4: Circulation decay and best-fitted functions at four different radial distances for 77 = 0.0702: solid 
line indicates e~ aT , where a is functions of rj and radial distance. 
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Figure 6 : Same as Fig. 4 but for 77 = 0.5031: dashed line indicates e bT , where b is functions of 77 and radial 
distance. 
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a Gaussian decay for weak and moderate turbulence 
but an exponential decay for strong turbulence. The 
model coefficients do not appear to much fluctuate 
with varying 77 in spite of the highly irregular na- 
ture of turbulence. In order to obtain a representa- 
tive model coefficient, more experiments have been 
performed for different 77 values and for different 
initial turbulence fields. The resulting coefficients 
(Fig. 8 ) fluctuate with varying 77 within the ranges 
of ci = 0.06 — 0.09 and c 2 = 0.08 - 0.20, respec- 
tively. The mean values of the coefficients, which 
can be considered as representative coefficients, are 
0.08 for model E and 0.13 for model G, respectively. 

Using the mean values of the coefficients, the 
model predictions of the average circulations for var- 
ious ambient turbulence levels are plotted in Fig. 9. 
The vortex lifespans in Fig. 9, defined as the time at 
which linking of a vortex pair occurs, are obtained 
from a model recently developed by Sarpkaya . 10 This 
model has been shown to agree well with data from 
water tank experiments 13 and LES predictions . 18 
Figure 9 indicates that until a linking of the vor- 
tices occurs the circulation is reduced by only about 
2 % of its initial value for weak ambient turbulence 
(77 = 0.03), while it may be reduced by as much as 
20 % for strong ambient turbulence (77 = 0.50). 



Figure 8 : Model coefficients with varying 77 for the 
circulation decay in Eg. (13) and Eg. (14) as obtained 
from numerical simulations. Different symbols rep- 
resent the results from different initial turbulence 
fields. Solid lines denote the average of the coeffi- 
cients for each model distinguished by dashed line. 


3.2 Vortex descent 


The descent of an ideal vortex pair, of which vor- 



Figure 9: Model prediction with varying 77 for the 
average circulation with vortex lifespan presented. 


ticity is concentrated in two parallel lines, will follow 
the relation of H = T in an inviscid fluid. In the nu- 
merical simulation, however, descent rate of a pair 
of the vortices can deviate from H/T = 1 due to 
the influence of image vortices implicit upon the as- 
sumption of periodic boundary conditions and the 
use of a vortex model different from ideal line vor- 
tices. From a two-dimensional laminar simulation 
with the same domain size as the crossplane size of 
the previous three-dimensional simulation and with 
a very small constant viscosity ( R e = T/u ~ 10 7 ), 
the descent rate intrinsic in the present domain size 
and vortex model appears to be H/T ~ 0.98, slightly 
less than ideal one. 

In order to develop a model for vortex descent 
due to ambient turbulence, we first assume that the 
vortex descent rate is proportional to the circulation 
of the accompanying vortex at r = bo, i.e., 


dh Ti 
dt 27T&0 ’ 


(17) 


where Ti is the circulation at r = bo. In the dimen- 
sionless form 


dH_ r, 
dT ~ Too 


(18) 


As discussed in the previous section, the circulation 
at r = b 0 is expected to decay as a Gaussian func- 
tion. Thus, substituting Eg. (16)** into Eg. (18), we 
obtain 


— ^ e ~(divT) 2 

dT 


(19) 


**At r = &o> £< 7 . (16) can be written as T — Toe c 2 T ? 2 ' r ’ 2 . 
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where d\ is an empirical coefficient proportional to 
c 2 and an approximation Ti(T = 0 ) « F^ is used. 
After integration, we finally obtain 

0 R7 

H= ™erf(d l r ] n (20) 

d\Tj 

where erf denotes an error function and the constant 
0.87 is an optimal value determined from our LES 
results. 

The above model for vortex descent has been ap- 
plied to our LES results for weak and moderate tur- 
bulence. As shown in Figs. lOa-d, the model predic- 
tions agree well with the LES data, although the co- 
efficient di appears to fluctuate with varying 77 . For 
strong turbulence, however, an optimal value for d\ 
has not been found. This may be due to the dom- 
inance of the exponential decay of the circulation 
over radial distances less than about 0.6 6 0 , imply- 
ing that the descent rate is not determined only by 
the circulations near r = 6 0 , but it is significantly in- 
fluenced by those at small radial distances less than 
6 0 • Nevertheless, for strong turbulence Eq.( 20) can 
still be matched with the LES results by fixing d\ 
and by changing the constant 0.87 to smaller values, 
as shown in Fig. lOe and f. For this case, therefore, 
the descent distance based on the LES results can 
be expressed by 

H = ollr r/<0 ■ 28 ’’ ^, ' <21) 

where d 2 is an empirical coefficient less than 0.87. 

The model coefficients di and d 2 for more cases 
have been obtained and plotted together in Fig. 11 
as done before for c\ and c 2 . As seen in Fig. 11 , the 
coefficients fluctuate within only maximum 17% of 
their mean values of 0.71 for model E and 0.84 for 
model G, respectively. 

Using the mean values of the coefficients, the 
model predictions for the vortex descent for various 
ambient turbulence levels are plotted in Fig. 12 again 
with the vortex lifespans obtained from Sarpkaya’s 10 
model. Figure 12 shows clearly that the descent 
rate of the vortices decreases with increasing ambi- 
ent turbulence due to the increasing rate of the cir- 
culation decay with increasing ambient turbulence. 
On the other hand, for weak turbulence, the de- 
scent distance appears to almost follow the ideal line 
H = 0.98 T before the vortex linking occurs. 

4. Comparison of observation data 
and model predictions 


In this section predictions of wake vortex decay 
and descent from the models proposed in the pre- 
vious section are compared to observation from the 
Memphis Field Program . 35 As part of the NASA 
sponsored field program, MIT Lincoln laboratory 
operated 10.6 fim C0 2 continuous wave laser Doppler 
radar (Lidar) and measured the line-of-sight veloc- 
ities of the vortices in a plane perpendicular to the 
flight path of approaching or departing aircraft. Cir- 
culation values were estimated from the measured 
line-of-sight velocity field . 35 

In this study, six cases (Table 2 ) are chosen based 
on atmospheric stability conditions and number of 
data points. Some of the above cases have been also 
analyzed by Sarpkaya 10 for a different purpose. The 
proposed models only address the effect of ambient 
turbulence on vortex decay and did not include ther- 
mal stratification. The initial conditions for the pro- 
posed models used the representative vortex separa- 
tion and far-field circulation in Table 2, as well as 
the TKE dissipation rate, e, measured at 40 m (Ta- 
ble 3). As seen from Table 3, the 77 values calculated 
by the definition (Eq.( 2)) are within the range of 
those used in our previous LES. 

Shown in Fig. 13 are the measured circulations 
averaged over a range of radial distances from 0.4 bo 
to 0.6 6 0 with model predictions based on Eq.(lS) 
and Eq.( 14) using the mean coefficients. The ini- 
tial average circulation ( Tq ) is determined from the 
present vortex model (Eq.(8)) with the theoretical 
circulation at r >> r c (Poo) given in Table 2. The 
data for times less than T — 1 has been discarded 
from the data set for all the cases analyzed, since 
the observed circulations in these times are subject 
to trigonometry errors (e.g., see Campbell et al. 35 ). 

As seen in Fig. 13, except for the cases in the sta- 
ble stratification (M-1252 and M-1273), predictions 
by the proposed models appear to be in reasonable 
agreement with observations, although data scatter 
is large. Significant differences in the rate of decay 
for the weak turbulence cases of M-1252 and M-1273 
between model predictions and observations may be 
caused by the interaction between vortex and sta- 
ble stratification, which is not included in numerical 
simulations. 

Also shown in Fig. 14 are the measured vortex 
descent with model predictions based on Eq. (20) 
and Eq. (21) again using the mean coefficients. In 
the calculation, the starting altitude of the vortices 
has been inferred from the Lidar measurement as 
in Robins et al ., 36 since the vortices appears to be 
located above the expected position, i.e., the gen- 


13 

American Institute of Aeronautics and Astronautics 


0 1 2 3 4 5 6 7 

T 


limensional simulations compared with the theoretical 


4 

onautics and Astronautics 









Flight 

number 

Aircraft 

type 

bo 

(m) 

Too 

Aircraft 
altitude (m) 

Atmospheric stability 
(observing time) 

M-1252 

B-757 

29.8 

294 

160.2 

stable (midnight) 

M-1273 

DC-10 

39.6 

400 

149.9 

stable (midnight) 

M-1409 

B-727 

25.8 

240 

146.6" 

neutral (evening) 

M-1569 

DC-9 

22.4 

215 

127.5 

neutral (evening) 

M-1581 

B-757 

29.8 

270 

166.1 

unstable (noon) 

M-1584 

DC-9 

22.4 

199 

125.7 

unstable (noon) 


Table 2: Initial parameters for wake vortices for the flights analyzed. 



u 


Figure 11: Same as Fig. 8 but for the vortex descent 
in Eq.( 20) and Eq.( 21). 



T 


Figure 12: Same as in Fig. 9 but for the vortex de- 
scent. 


Flight 

number 

€ 

(m 2 s -3 ) 

1 

M-1252 

0.55 x 10" 5 

0.03 

M-1273 

0.20 x 10~ 3 

0.12 

M-1409 

0.29 x lO - * 

0.28 

M-1569 

0.79 x 10~ 3 

0.17 

M-1581 

0.13 x 10" 1 

0.51 

M-1584 

0.10 x 10" 1 

0.43 


Table 3: Turbulence energy dissipation rate mea- 
sured at elevation of 40 m and corresponding dimen- 
sionless turbulence intensity. 

erating aircraft height in Table 2. Consistent with 
the large underestimation for the circulation decay, 
the model predictions in the cases of M-1252 and 
M-1273 significantly overestimate the observed vor- 
tex descent especially at later time period. Except 
for the case of M-1409 (Fig. 14c), the model predic- 
tions for the other cases appear to be in reasonable 
agreement with observations. The underprediction 
in the case of M-1409 could be caused by other fac- 
tors such as wind shear or large-scale atmospheric 
vertical motion. 

5, Summary and conclusions 

This study represents a first step at understand- 
ing how wake vortices behave in the atmosphere by 
isolating the effects of ambient turbulence on vor- 
tex decay. LES simulations are conducted by first 
growing a field of ambient turbulence which is nearly 
isotropic and spatially homogeneous. A vortex pair, 
representative of aircraft trailing vortices, is intro- 
duced once the turbulence field reaches a steady 
state. Results show that the decay rate of the vortex 
circulation increases clearly with increasing levels of 
ambient turbulence. In these simulations Crow link- 
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Figure 13: Comparison of field observations and model predictions for the average circulation; (a) and (b) 
stable environment; (c) and (d): neutral environment; (e) and (f): unstable environment. 
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ing and its subsequent effect on vortex decay are 
suppressed by the choice of domain size. 

Based on the LES results, simple decay mod- 
els are proposed as functions of dimensionless ambi- 
ent turbulence intensity (77) and dimensionless time 
(T) for the circulation averaged over radial distances 
from 0.4 60 to 0.6 &o- With good agreement with the 
LES data, a Gaussian type of vortex decay model is 
applicable for weak and moderate turbulence, while 
an exponential type of vortex decay model can be 
applied for strong turbulence. The LES simulations 
and the fitted models show that the turbulence dissi- 
pation rate characterizes the level of ambient turbu- 
lence responsible for vortex decay. A model for the 
vortex descent based on above vortex decay model 
is also proposed as functions of 77 and T and can 
be represented by an error function. The proposed 
models for the vortex decay and descent are applied 
to available field data obtained from the Memphis 
airport. In a turbulent atmosphere with either neu- 
tral or unstable stratification, the model predictions 
appear to be in reasonable agreement with the ob- 
servation data. In a stably-stratified atmosphere 
with low turbulence, the models largely underesti- 
mate the observed circulation decay with consistent 
overestimation of the observed vortex descent. The 
lack of agreement in the latter is attributed to the 
effects of stratification on vortex decay. 
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Appendix A 


According to Donaldson and Bilanin , 7 hereafter 
referred as DB, for the neutrally stratified and high 
Reynolds number flows, the equation that governs 
the rate of change for mean axial vorticity £ due to 
turbulent motion is given by 


dC 

dt 




Ml) 


transport assumption, 


«C) = (A5) 

where v e is an effective eddy viscosity. The radial 
gradient of mean axial vorticity may be approxi- 


mated by 

dr r 

(46) 

and 

r 



C ~ ~2 > 

(47) 

so that 




d c ^ r 

dr r 3 ’ 

(48) 


Note that dQ /dr is assumed to be proportional to 
the difference of mean axial vorticity at the vortex 
center and that at radial distance r. Thus, Eq.(A4) 
can be expressed as 


dT v e Y 

dt r 2 


(A9) 


Above derivations are almost the same as DB’s ex- 
cept for the fact that the dependency of circulation 
decay on the radial distance r is considered in the 
present study, while bo is used instead of r in DB’s 
study, since DB consider the decay rate of the cir- 
culation over an oval within which a pair of vortices 
resides. 

From an equilibrium second-order closure theory, 
DB obtain the relation 


Integration of Eq.(Al) over region defined by the 
radial distance r yields 

SH dA= -H^ iA - w 

The area integral on the right-hand side of Eq.(A2) 
can be transformed into a line integral taken around 
the perimeter of the area using the Stoke’s theorem, 
i.e., 

^ ~ u ' n C' ds, (43) 

where u' n is the fluctuation of velocity normal to the 
surface. The line integral leads to 

^ = — 2 tt r(u^C), (44) 

where ( ) denotes an average around the perime- 
ter with the radial distance r. Under down-gradient 


v e = 0.26A<?, (A 10 ) 

where A is a measure of the integral length scale of 
turbulence. Assuming that 

A~|, (411) 

DB then obtain £<7.(1). However, their approxima- 
tion of £< 7 . (All) is not physically correct, since 60 
has no relation with background atmospheric tur- 
bulence whose length scale varies significantly over 
day and night. As already discussed in the intro- 
duction, in this study we use e, rather than q and 
A, as the ambient turbulence parameter relevant to 
circulation decay. 

For strong turbulence but at smaller radial dis- 
tances where stronger radial gradient of mean axial 
vorticity may exist, the parameterization of £< 7 -(A 5 ) 
would be acceptable. The effective eddy viscosity in 


20 

American Institute of Aeronautics and Astronautics 


this case can be expressed in terms of effective tur- 
bulence velocity and length scales, i.e., 

»e ~{eb 0 ) 1/3 b 0 . (A12) 


Note that (ebo) 1 ^ 3 is a characteristic turbulence ve- 
locity scale at the scale of the vortex separation dis- 
tance b 0 . Substituting Eq.(Al2) into Eq.(A9), we 
obtain 


dr 

dt 


(eb 0 ) 1/3 b 0 


r. 


(> 113 ) 


Thus, the decay rate of the circulation depends on 
both the turbulence energy dissipation rate e and ra- 
dial distance r, unlike DB’s model (Eq.( 1)). Equa- 
tion (A13) can be nondimensionalized as 




dT ~ R 2 * 


(A14) 


For weak and moderate turbulence or at larger 
radial distances for strong turbulence where the ra- 
dial gradient of mean axial vorticity may be weak, 
the parameterization of Eq.(A5) would not be ac- 
ceptable. It is anticipated that the radial turbu- 
lent vorticity flux u ' n £' in this case will be mostly 
influenced by turbulent diffusion of vorticity from 
the core region rather than by local diffusion like in 
Eq.(A5). This implies that would be time de- 
pendent. We can hypothesize that the change rate 
of in this case depends on turbulence velocity 
and mean axial vorticity gradient from vortex center 
to given r, i.e., 


d«C') 

dt 


“/((‘M'".?)- (-U5) 


(AW) 


From dimensional arguments, 

aKC 7 ) W 2/3 £ 

dt r 3 

Integration of Eq.(A16) with respect to time leads 
to 

(e&o) 2/3 <r 


<c ~ -■ 


(> 117 ) 


where the initial value of u' n (,' is assumed to be zero. 
Substituting Eq.(A17) into Eq. (.44), we obtain 


dT 

dt 


(eb 0 ) 2/3 t 


r. 


(>118) 


The above equation can be nondimensionalized as 
dT q 2 T 


dT 


R 2 


r. 


(A19) 
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